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Abstract. The adhesion zone of immune cells, the 'immunological synapse', exhibits 



characteristic domains of receptor-ligand complexes. The domain formation is 

likely caused by a length difference of the receptor-ligand complexes, and has been 

C^ investigated in experiments in which T cells adhere to supported membranes with 

anchored ligands. For supported membranes with two types of anchored ligands, 
MHCp and ICAM1, that bind to the receptors TCR and LFA1 in the cell membrane, 
the coexistence of domains of TCR-MHCp and LFA1-ICAM1 complexes in the cell 
adhesion zone has been observed for a wide range of ligand concentrations and affinities. 
For supported membranes with long and short ligands that bind to the same cell 
receptor CD2, in contrast, domain coexistence has been observed for a rather narrow 
ratio of ligand concentrations. In this article, we determine detailed phase diagrams 
for cells adhering to supported membranes with a statistical-physical model of cell 

("**•) adhesion. We find a characteristic difference between the adhesion scenarios in which 

two types of ligands in a supported membrane bind (i) to the same cell receptor or (ii) 
to two different cell receptors, which helps to explain the experimental observations. 

C ■ Our phase diagrams fully include thermal shape fluctuations of the cell membranes 

on nanometer scales, which lead to a critical point for the domain formation and to a 



cooperative binding of the receptors and ligands. 
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1. Introduction 

Cell adhesion is mediated by the specific binding of a variety of membrane- anchored 
receptor and ligand molecules. In 1990, Springer suggested that the length difference of 
receptor-ligand complexes in the contact zone of immune cells may lead to segregation, 
i.e. to the formation of domains within the cell contact zone that contain receptor- 
ligand complexes with different lengths fl]. The 'length' of a receptor-complex here is 
the intermembrane distance, or local membrane separation at the site of the complex. 
A length difference between receptor-ligand complexes leads to an indirect, membrane- 
mediated repulsion of the complexes because the membranes have to bend to compensate 
the mismatch, which costs bending energy. Important receptor-ligand complexes in T- 
cell adhesion are the TCR-MHCp complex with a length of about 13 nm p], the CD2- 
CD48 complex with the same length of 13 nm [3}j5], and the LFA1-ICAM1 complex with 
a length of about 40 nm J6l. In 1998 and 1999, the contact zone of T cells was indeed 
found to contain domains that either contain the short TCR-MHCp or the long LFAl- 
ICAM1 complexes [7j[8]. As expected from their length, the CD2-CD48 complexes 
are located within the TCR-MHCp domains |5]. However, the question whether the 
domain formation is predominantly caused by the length mismatch of receptor-ligand 
complexes is complicated by the role of the actin cytoskeleton, which polarizes during 
T-cell adhesion and transports clusters of TCR-MHCp complexes towards the center of 
the cell contact zone (9 11, and by additional, direct protein-protein interactions 12 



The domain formation is closely linked to T-cell activation, with TRC clusters forming 
within seconds of T-cell adhesion triggering the first activation signals |9 13 . 



Direct evidence for a central role of the length of receptor-ligand complexes comes 
from experiments in which these lengths are altered by protein engineering (5, 14 



Milstein and coworkers [5] have considered variants of the protein CD48 with four and 
five immunoglobolin-like (Ig-like) domains. The CD48 variants are longer than the 
CD48 wildtype, which contains only two Ig-like domains. The CD48 wildtype and both 
CD48 variants bind to CD2 on T cells. From electron micrographs of the contact zone 
between T cells and supported membranes that contain one of the three CD48 types, 
Milstein and coworkers found that the length of the CD2-CD48 complex is 12.8 ± 1.4 
nm for wildtype CD48, 14.2 ± 1.2 nm for the CD48 variant with four Ig-like domains, 
and 15.6 ± 1.4 nm for the variant with five Ig-like domains. In fluorescence experiments 
of T cells on supported membranes that contain mixtures of two of the three CD48 
types, Milstein and coworkers observed that CD2-CD48 wildtype complexes segregate 
from both CD2-CD48 variant complexes. The segregation seems to be driven by the 
length difference of the complexes since the T-cell cytoskeleton can only 'act on' CD2 
and, thus, can hardly 'discriminate' between the different complexes. However, Milstein 
and coworkers observe domain coexistence in the contact zone only within a narrow 
range of concentration ratios of CD48 wildtype and CD48 variants. For T cells adhering 
to supported membranes with MHCp and ICAM1, in contrast, domain coexistence 
has been observed for a rather wide range of MHCp and ICAM1 concentrations and 
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In this article, we calculate detailed phase diagrams for cells adhering to supported 
membranes with anchored ligands. We consider two general adhesion scenarios: In the 
first scenario, long and short ligands in the supported membrane bind to the same cell 
receptor (see section 4), as in the experiments of Milstein and coworkers (5J, in which 
CD48 wildtype and a CD28 variant in the supported membrane both bind to CD2 in 
the T-cell membrane. In the second scenario, two types of ligands in the supported 
membrane bind to two types of receptors in the cell membrane (see section 5), as in 
experiments in which MHCp and ICAM1 in the supported membrane bind to TCR and 
LFAl in the T-cell membrane. We find a characteristic difference between the phase 
diagrams in the two scenarios (see fig. [5]). In the first scenario, domain coexistence only 
occurs along a coexistence line. In the second scenario, in contrast, domain coexistence 
occurs in a wide coexistence region. Our phase diagrams thus help to understand 
why Milstein and coworkers observe domain coexistence only within a narrow range 
of concentration ratios. 

Our calculations are based on a statistical-physical model of cell adhesion (see 
section 2). In this model, the membranes are described as elastic sheets discretized into 



small patches that can contain single receptor or ligand molecules 16 19 . The binding 
and domain formation of receptor-ligand complexes is affected by thermally excited 
shape fluctuations of the membranes on nanometer scales. These shape fluctuations 
lead to a critical point for the segregation of long and short receptor-ligand complexes 



19 20 . The critical point depends on the length difference of the complexes, on the 
concentrations and affinities of the receptors and ligands, and on the bending rigidity of 
the membranes (see section 4). The critical point constitutes a threshold for segregation, 
or domain formation, and may help to understand why Milstein and coworkers have 
observed segregation of wildtype CD48 from each of the two CD48 variants, but not 
segregation of the two CD48 variants [5]. In addition, the membrane shape fluctuations 
on nanoscales lead to a cooperative binding of receptor-ligand complexes [21] (see section 
3). 

2. Statistical-physical description of cell adhesion 

Cell adhesion involves length scales that differ by orders of magnitude (see fig. [I]) . The 
diameters of the cell and cell contact zone have values of several micrometers, while 
the average separation of the membranes within the contact zone is typically tens of 
nanometers. Other important length scales in the cell contact zone are the average 
distance between receptor-ligand bonds, and the binding width of receptor and ligand 
molecules. The binding width is the difference between the smallest and the largest 
local membrane separation at which the molecules can bind. The binding width of the 
typically rather stiff receptor and ligand proteins that mediate cell adhesion is much 
smaller than the length of the proteins. 

The binding equilibrium and segregation of receptor-ligand complexes in cell contact 
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zones is affected by membrane shape deformations and fluctuations. Since bound 
receptor-ligand complexes constrain the local separation of the membranes, the relevant 
deformations and fluctuations of the membranes occur on lateral length scales smaller 
than the average distance between neighboring pairs of complexes, which is about 100 
nanometers for complex concentrations of about 100 per square micrometer J8j. It 
is reasonable to assume that the elasticity of the membranes is dominated by their 
bending rigidity on these length scales. The binding rigidity k dominates over the 
membrane tension a on lateral length scales smaller than the crossover length \/n/a 



22 , which is of the order of several hundred nanometers for cell membranes 23 



The cytoskeletal elasticity 24-27 contributes on length scales larger than the average 



distance between the cytoskeletal anchors in the membrane, which may be around 100 



nanometers 28 . The bending rigidity thus is likely to dominate over the lateral tension 
and the cytoskeletal elasticity on lateral length scales up to 100 nanometers relevant 
here. 

We have developed discrete models for the adhesion of membranes via anchored 



receptors and ligands 16,18,19,29 . In discrete models, the two apposing membranes in 



the contact zone of cells or vesicles are divided into small patches 16 18 20 23 29 36 



In our models, the rigidity-dominated elasticity of the membranes in the contact zones 



of cells or vesicles is described by [16, 18 

^el{0 = ^ E ( A <^) 2 C 1 ) 

i 

where U is the local separation of the apposing membrane patches i. The elastic energy 
depends on the mean curvature (Adh)/a 2 of the separation field U with the discretized 
Laplacian Adk = ^1+^2+^3+^4- 4/j. Here In to la are the membrane separations at the 
four nearest-neighbor patches of membrane patch i on the quadratic array of patches. 
The linear size a of the membrane patches is chosen to be around 5 nm to capture the 
whole spectrum of bending deformations of the lipid membranes (37|. The 'effective 
bending rigidity' of the two membranes with rigidities K\ and Ki is k = ^1^2/(^1 + ^2)- 
If one of the membranes, e.g. membrane 2, is a planar supported membrane, the effective 
bending rigidity k equals the rigidity K\ of the apposing membrane since the rigidity k 2 
of the supported membrane is taken to be much larger than K\. 
The overall energy of the membranes in the cell contact zone 

1-L{l,n,m} =H e \{l} + H mt {l,n,m} (2) 

is the sum of the elastic energy % e i{l} and interaction energy 7iwt{l,n,m}. The 
interaction energy depends on the distribution n of the receptors in membrane 1, on 
the distribution m of the receptors in membrane 2, and on the separation field / of the 
membranes. In our models, each patch of the discrete membranes can only be occupied 
by one receptor or ligand molecule. Mobile receptor and ligand molecules diffuse by 
'hopping' from patch to patch, and the thermal fluctuations of the membranes are 
reflected in variations of the local separation of apposing membrane patches. A receptor 
can bind to a ligand molecule if the ligand is located in the membrane patch apposing 
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Figure 1. A cell adhering to a supported membrane with anchored ligands that 
bind to receptors in the cell membrane. The binding of receptors and ligands in the 
cell adhesion zone is affected by membrane shape deformations and fluctuations on 
nanometer scales, which are dominated by the bending rigidity of the cell membrane. 
The immune cell receptors are typically mobile along the membrane and not, or only 



weakly 11 , coupled to the cytoskeleton. 



the receptor, and if the local separation of the membranes is close to the length of the 
receptor-ligand complex (see fig. [2] and below). In discrete models, the receptor and 
ligand molecules are taken into account as individual molecules. In continuum models, 
in contrast, the distributions of receptor and ligand molecules on the membranes are 



described by continuous concentration profiles 38-48 



3. Adhesion via a single type of receptor-ligand complexes 

3.1. Interaction energy of receptors and ligands 

We first consider the case in which the adhesion is mediated by a single type of receptor- 
ligand complexes. Examples of this case are (i) cells adhering to supported membranes 
that contain a single type of ligand (5, 49-51], and (ii) vesicles with anchored receptors 



that adhere to supported membranes or surfaces with complementary ligands 35 ,52 62 . 



Segregation of receptor-ligand complexes in cell adhesion zones 




fefe 






I 

I 




/» 



Figure 2. A supported membrane with anchored ligands (bottom) that bind to 
receptors in an apposing cell or vesicle membrane (top). In our model, the membranes 
are discretized in small patches, which can contain single receptor or ligand molecules. 
The shape and thermal fluctuations of the cell or vesicle membrane in the adhesion 
zone are described by the local separations li of apposing membrane patches i. A 
receptor can bind to a ligand molecule (i) if the ligand is located in the membrane 
patch apposing the receptor, and (ii) if the local membrane separation li is close to the 
length of the receptor-ligand complex. 



The interactions of receptors and ligands within the contact zone of the cell or vesicle 
are described by the interaction energy [17,21 

Wint{l, n, m} = 22 nimiV{li) (3) 

i 

in our model. Here, the occupation number n t = 1 or indicates whether a receptor is 
present or absent in membrane patch i of the cell, and rrii = 1 or indicates whether a 
ligand is present or absent in patch % of the apposing membrane. Receptor and ligand 
molecules in apposing patches % of the membranes interact with the potential V(lj). For 
simplicity, we describe this interaction by the square-well potential 

V{li) = - U for l - Z we /2 < k < l + l we /2 

= otherwise (4) 

which depends on the binding energy U > 0, and the length l and binding width Z we 
of a receptor-ligand complex. A receptor thus binds to an apposing ligand with energy 
—U if the local separation 1$ of the membranes is within the binding range l ± Z we /2. 



3.2. Effective adhesion potential 

The binding equilibrium of the membranes in the contact zone can be determined from 
the free energy J 7 = —ksTlnZ, where Z is the partition function of the system, ks is 
Boltzmann's constant, and T is the temperature. The partition function Z is the sum 
over all possible membrane configurations, with each configuration {I, n, m} weighted 
by the Boltzmann factor exp [— 7{{l,n, m}/fcsT]. A membrane configuration in the 
contact zone is specified by the separation field / of the membranes, the distribution 
m of the receptors in the cell membrane, and the distribution n of ligands in the 
apposing membrane. In our model, the partial summation in the partition function 



Segregation of receptor-ligand complexes in cell adhesion zones 7 

Z over all possible distributions m and n of receptors and ligands can be performed 
exactly, which leads to an effective adhesion potential. The effective adhesion potential 
V e f(li) is a square- well potential with (i) the same binding range l we as the receptor- 
ligand interaction ^ and (ii) an effective potential depth U e f that depends on the 



concentrations and binding energy U of receptors and ligands [17 18,21 : 

V ei (k) = - U ei for l - l wc /2 <k<l + l wc /2 

= otherwise (5) 

For typical concentrations of receptors and ligands in cell adhesion zones up to hundred 
or several hundred molecules per square micrometer, the average distance between 
neighboring pairs of receptor and ligand molecules is much smaller than the width 
of the molecules. For these small concentrations, the effective binding energy of the 



membranes is 21 



U e{ nk B Ta 2 e u/kBT [R][L] (6) 

where [R] is the area concentration of unbound receptors in the cell membrane, and [L] 
is the area concentration of unbound ligands in the apposing membrane. The binding 
equilibrium in the contact zone thus can be determined from considering two membranes 
with the elastic energy (fiT) that interact via an effective adhesion potential with well 
depth U e f and width Z wc . 

3. 3. Area fraction P& of the membranes within binding range of receptors and ligands 

Receptor-ligand complexes can only form at membrane patches with a local separation 
within the binding range l ± l wc /2 of the receptors and ligands (see eq. (J5l)). The area 
concentration [RL] of the receptor-ligand complexes in the contact zone therefore is 



proportional to the fraction p, of these membrane patches 21 



[RL]*P b K[R][L] (7) 

Here, K is the equilibrium constant for receptor-ligand binding within this membrane 
fraction. In our model, the equilibrium constant is K = a 2 e u ^ kBT . 

In equilibrium, the fraction P& of membrane patches with a local separation within 
receptor-ligand binding range depends on the effective binding energy U e f, the binding 
width Z we , the effective rigidity k of the membranes, and the temperature T. We 
have found that the effect of these four quantities on p, can be captured by a single 



dimensionless quantity, the rescaled effective potential depth 21 

u = U e{ K llJ{k B Tf « {K/k B T)ll c K[R} [L] (8) 

To a first approximation, the membrane fraction p, depends only on u for typical 
lengths and concentrations of receptor-ligand complexes in cell adhesion zones. In 
cell adhesion zones, direct contacts between the membranes can be neglected since the 
average separation of the membranes, which depends on the length l of the complexes, 



is typically larger than the thermal membrane roughness 21 . 
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Figure 3. An important quantity is the area fraction Pf, of the membranes within 
binding separation of the receptors and ligands. The area fraction Pi, (shown in 
red) increases with the concentrations of receptors and ligands, since the formation 
of receptor-ligand bonds 'smoothens out' thermal membrane shape fluctuations. The 
'smoothening' facilitates the formation of additional receptor-ligand bonds and, thus, 



leads to a binding cooperativity 21 



From Monte Carlo simulations, we have found that the functional dependence of 
the area fraction P^ on the rescaled potential depth u is well described by 

P ~ u 
f b ■ — 

Ci+U 

with the dimensionless coefficient c\ 



0.071 



21 



(9) 



The membrane fraction Pu increases 



with u and, thus, increases with the effective binding energy U e { and the effective bending 
rigidity k. The reason for this increase is that the roughness of the membranes resulting 
from thermal shape fluctuations decreases with C/ e f and k. The membrane fraction P& 
decreases with the temperature T since the roughness increases with T. The thermal 
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roughness, defined as the standard deviation of the local membrane separation from its 
average, is the characteristic length scale for membrane excursions in the perpendicular 
direction. The membrane fraction P& within receptor-ligand binding range is much 
smaller than 1 if the roughness is large compared to the binding width Z wc of the 
complexes, and close to 1 if the roughness is small compared to l we . 

3.4- Concentrations of bound and unbound receptors of an adhering cell 
From eqs. (J7|) to (J9l), we obtain the relation 21 



lRL] „ <jnmi no) 

[RLl ~ Cl k B T + nll c K[R][L] (10) 

between the area concentration [RL] of bound receptor-ligand complexes in the contact 
zone and the area concentrations [R] and [L] of unbound receptors and ligands. This 
nonlinear relation reflects the cooperative binding of receptors and ligands. This 
cooperativity arises because the binding of receptors and ligands suppresses thermal 
membrane fluctuations and, thus, smoothens the membranes, which facilitates the 
binding of additional receptors and ligands (see fig. |3j. 

The total number N of receptors in the cell membrane is constant. The 
concentrations of bound and unbound receptors are therefore connected by the 
additional relation 119 



N&[R]A+[RL]A C (11) 

where A is the total area of the cell membrane, and A c the contact area. We have 
neglected here the area occupied by bound receptor-ligand complexes since this area is 
small compared to the total contact area A c for typical concentrations in cell adhesion 
zones. The concentrations of unbound receptors within and outside of the contact area 



then are equal. Together, the two relations (10) and ([111) determine the concentration 
[R] of unbound receptors and the concentration [RL] of bound receptors in the contact 
zone. 

4. Two types of membrane-anchored ligands adhering to the same cell 
receptor 

4-1. Interaction energy of receptors and ligands 

In recent experiments by Milstein and coworkers J5], long and short ligands anchored to a 
supported membrane bind to the same receptor of an adhering T cell. These ligands are 
wildtype CD48 and elongated CD48 variants, and the receptor in the T cell membrane 
is CD2. In our model, this situation is described by the interaction energy 

n int {l,n,m} = Y^MSmtMh) + S mi>2 V 2 (h)) (12) 

i 

for the two apposing membranes in the cell contact zone. Here, the occupation number 
rrii = 1, 2, or indicates whether a ligand L 1 of type 1, a ligand L 2 of type 2, or 
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no ligand is present in patch i of the supported membrane, and Ui = 1 or indicates 
whether a receptor R is present or not in the apposing patch i of the cell membrane. 
The Kronecker symbol 5ij equals 1 for i = j and is equal to for % ^ j. The potential 
V\ thus describes the interaction of the receptor R with the ligand protein L\ , and the 
potential V 2 the interaction between R and L 2 . For simplicity, V\ and V 2 are again taken 
to be 

Vi(U) = U x for h - / we /2 < k < h + l we /2 

= otherwise (13) 

and 

V 2 (k) = U 2 for l 2 - / wc /2 < k <l 2 + / wc /2 

= otherwise (14) 

with binding energies U\ and U 2 and equilibrium lengths h < l 2 of the complexes RL\ 
and RL 2 . We have assumed here that the two complexes have the same binding width 

'we' 

4-2. Effective adhesion potential 



As in section |3.2[ the summations over all possible distributions m and n of receptors 

and ligands in the partition function of the model leads to an effective adhesion potential 

19 20 . The effective adhesion potential now is a double- well potential (see fig. [4]). Both 



wells have the same width Z we as the potentials (13) and (14). The well with its center 
at the membrane separation U = l\ reflects the interactions of the receptors R with 
the shorter ligands L l5 and the well centered at U = l 2 reflects the interactions of the 
receptors and the longer ligands L 2 . In analogy to eq. (|6]), the depth of the two wells 

Uf ^ k B T KxiRJiLx] (15) 

and 

U^k B TK 2 [R][L 2 ] (16) 

depend on the concentrations [R], [L\] and [L 2 ] of unbound receptors and ligands, and 
on the binding constants K\ = a 2 e Ul ^ kBT and K 2 = a 2 e U2 ^ kBT for receptors and ligands 



within the appropriate binding ranges 19 , 20 . The binding equilibrium of the receptors 



R and ligands L\ and L 2 in the contact zone thus can be determined from considering 
two apposing membranes with elastic energy rtl]) that interact via an effective double-well 



potential with well depths U± and U 2 given by eqs. (15) and (16). 



4-3. Phase diagram 

If the two wells of the effective adhesion potential are relatively shallow, thermal 
membrane fluctuations can easily drive membrane segments to cross from one well to the 
other. If the two wells are deep, the crossing of membrane segments from one well to the 
other well is hindered by the potential barrier of width ly, a between the wells (see fig. [4]). 
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Figure 4. (a) A supported membrane with short (green) ligands L\ and long (red) 
ligands L-i that bind to the same receptor R in the cell membrane (top). - (b) The 
interactions of the receptors and ligands lead to an effective double-well adhesion 
potential V c { of the membranes. The potential well 1 at small membrane separations I 
reflects the interactions of the receptors with the short ligands, and the potential well 2 
at larger membrane separations the interactions with the long ligands. The depths Uf 
and C/f f of the two potential wells depend on the concentrations and binding constants 
of the receptors and ligands (see eqs. ( fl5| ) and (16])). The wells have the same width Z we 
as the receptor-ligand interactions ( |13[ ) and (14), and a separation Zba = h-~h+lwc that 
depends on the difference between the equilibrium lengths l\ and li of the complexes 
RL\ and RL2- 



The potential barrier induces a line tension between adjacent membrane segments that 
are bound in different wells 



63 



Beyond a critical depth of the potential wells, the 
line tension leads to the formation of large membrane domains that are bound in well 
one or well two. Within each domain, the adhesion of the membranes is predominantly 
mediated either by the receptor-ligand complexes RL\ or by the complexes RL2. 

We have previously found that the critical potential depth for domain formation is 

c(k B T) 2 



U. 



cf 



kJ 1 



(17) 



we'-ba 



with the prefactor c = 0.225±0.02 determined by Monte Carlo simulations 20 . Domain 



formation in the contact zone or, in other words, segregation of the complexes RL\ and 
RL 2 can only occur if the effective potential depths Uf and Uf exceed the critical 
potential depth Uf. The critical potential depth depends on the temperature T and 
the bending rigidity k as well as on the width Z we and separation Zba of the two potential 
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wells. In deriving eq. (17), we have neglected direct membrane-membrane contacts, 



which is reasonable for typical concentrations and lengths of receptor-ligand complexes 



in cell adhesion zones 20 21 . For these complex concentrations and lengths, the thermal 
membrane roughness is smaller than the lengths of the receptor-ligand complexes. 
Domain coexistence occurs for equal depths 

Uf = Uf (18) 



of the potential wells if the two wells have the same width l we as in Fig. 4 With eqs. ( 15 ) 



and (16), this coexistence condition implies that domain coexistence occurs along the 
line with 

K 1 [L 1 ] = K 2 [L 2 ] (19) 

in the [Li]-[L2J-plane. The line has the slope K\jK 2 and ends at a critical point (see 



phase diagram in fig. p[a)). For Uf > £/f f , we have ifi|Xi] > K 2 [L 2 ]. The adhesion 
is then dominated by the short complexes RL\ throughout the cell contact zone. For 
Uf < U 2 , in contrast, the adhesion is dominated by the long complexes RL 2 in the 
whole contact zone. If the supported membrane is much larger than the cell contact zone, 
it seems reasonable to assume that the concentrations [L\] and [L 2 ] of unbound ligands 
do not change upon adhesion since the 'ligand reservoir' in the supported membrane is 
large. The concentrations [L-\\ and [L 2 ] in our model then correspond to the experimental 
ligand concentrations in the supported membrane prior to adhesion, and our phase 
diagram in fig. J5](a) to the phase diagram in fig. 7 of Ref. (H], see Discussion. 

5. Two types of membrane-anchored ligands adhering to different cell 
receptors 

5.1. Interaction energy of receptors and ligands and effective adhesion potential 

Several experimental groups have investigated the adhesion of T cells to supported 
membranes with anchored MHCp and ICAM1 ligands |||9||ll||l5||6§]68|. The ligand 
MHCp binds to the T cell receptor (TCR), and the ligand ICAM1 to the integrin LFAl 
in the T cell membrane. The TCR-MHCp complex has a length of around 13 nm J2], 
and the LFA1-ICAM1 complex a length of 40 nm J6J. A situation in which two ligands in 
the supported membrane bind to different receptors in a cell membrane can be described 
in our model via the interaction energy [20] 



H int {l,n,m} = J2 (<W<W^i(*i) + S ni , 2 S mit2 V 2 (k)) (20) 

i 

Here, the occupation number n« = 1, 2, or indicates whether a receptor R\, a receptor 
R 2 , or no receptor is present in patch % of the cell membrane in the contact zone, while 
rrii — 1, 2, or indicates whether a ligand L±, a ligand L 2 , or no ligand is present in the 
apposing patch % of the supported membrane. The interaction of a receptor R\ with an 
apposing ligand L\ is described by the potential Vy(li), and the interaction of R 2 with L 2 



by the potential V 2 (k). As in section 4.2, a summation over all possible distributions n 
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Figure 5. (a) A supported membrane with two types of ligands, L\ and L2, that 
bind to the same receptor R of an adhering cell (top). The different length of the 
ligands causes a membrane-mediated repulsion between the receptor-ligand complexes 
RLi and RL2- For sufficiently large concentrations \L\\ and [L2] of the ligands, the 
repulsion of the complexes RL\ and RL2 leads to the formation of domains. However, 
domain coexistence in the cell adhesion zone only occurs for equal effective binding 
strengths Ki[Li] = K 2 [L2\ of the ligands. Here, K^ and K 2 are the binding equilibrium 
constants of the two ligands at appropriate membrane separations. Domain coexistence 
thus occurs along the shown line with slope Ki/K 2 in the [Li ]-[L 2 ] -plane (bottom). 
The line ends at the critical point for domain formation. For Jfj[ii] > .Kg [£2]) the 
adhesion is dominated by the ligand Li throughout the cell adhesion zone, and by the 
ligand L 2 for .Kx|Xi] < K 2 [i 2 ]. We have assumed here that the supported membrane 
is large compared to the adhesion zone, which implies that the concentrations [L\] 
and [L2] of unbound ligands do not change significantly upon adhesion. - (b) A 
supported membrane with two types of ligands, L\ and L2, that bind to different 
cell receptors i?i and i?2- Domain coexistence in the cell adhesion zone now occurs 
for ifi[i?i][Li] = K2 [R2] [L2] , which leads to a broad coexistence region in the [£1]- 
[L2]-plane since the concentrations of unbound receptors [Ri] and [R2] depend on the 
numbers of bound receptors and, thus, on the adhesion zone fractions occupied by the 
two domains. 



Segregation of receptor-ligand complexes in cell adhesion zones 14 

and m of receptors and ligands in the partition function leads to an effective double-well 
potential of the membranes. The effective potential has the same form as in fig. Kb), 
but the depths of the two wells 

Uf & k B T[Ri][Li]Ki (21) 

and 

Uf ^ k B T[R 2 ][L 2 ]K 2 (22) 

now depend on the concentrations [Ri] and [R 2 ] of unbound receptors in the cell 
membrane, on the concentrations [Li] and [L 2 ] of unbound ligands in the supported 
membrane, and on the binding constants K\ = a 2 e c/l//fesT and K 2 = a 2 e U2 ^ kBT for the 
complexes R\Li and R 2 L 2 |19|. 

5.2. Phase diagram 



As in section |4.3[ domain coexistence in the cell contact zone requires equal depths 
Uf = C/f of the potential wells if the two wells have the same width l we . The effective 
adhesion potential then is a symmetric double-well potential. We assume here again 
that the total area of the supported membrane is large compared to the cell contact 
zone. In this case, the numbers of bound ligands in the contact zone is negligible 
compared to the numbers of unbound ligands in the total supported membrane, which 
implies that the concentrations [L\\ and [L 2 ] of unbound ligands do not change during 
adhesion. However, the concentrations [Ri] and [R 2 ] of unbound receptors in general 
change during cell adhesion because the contact area is typically a substantial fraction 
of the overall area of the cell membrane, and because the total numbers N\ and N 2 of 
the receptors in the cell membrane are constant. During adhesion, a smaller or larger 
fraction of the receptors will form bound complexes R±Li or R 2 L 2 (see also section 



3.4). The concentrations [i?iLi] and [R 2 L 2 ] of the receptor-ligand complexes in the cell 
contact area depend on the fractions P\ and P 2 of the membranes within well 1 and 
well 2 of the effective adhesion potential. The receptor-ligand complexes R\L\ can only 
form in the membrane fraction Pi of the contact area within binding range l\ ± l we /2 of 
Ri and L\, and the complexes R\L\ only in the membrane fraction P 2 within binding 
range l 2 ± l wc /2 of R 2 and L 2 . 

For the symmetric double-well potential with Uf = Uf, the membrane fractions 
Pi and P 2 within well 1 and well 2 depend primarily on the rescaled potential depth 

u = UfKl 2 w J(k B T) 2 = UfnllJ{k B Tf (23) 

as in section 



3.3 The Monte Carlo data in fig. p[a) illustrate how Pi and P 2 depend 
on u. Below the critical potential depth u c , the membrane fluctuates between the two 
wells. Because of the symmetry of the potential, Pi and P 2 attain the same value 
Pb(u) = P\{u) = P 2 {u) for u < u c . Above the critical potential depth u c , we have a 
spontaneous symmetry breaking of the membranes into domains that are predominantly 
bound in well 1 or well 2, or in other words, predominantly bound by the complexes RiLi 
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or the complexes R2L2. The symmetry breaking for u > u c is reflected by two branches 
P b {u) and P b ~(u) of the membrane fraction P b (u) within the wells (see fig.pfa)). For the 
domain predominantly bound in well 1, we have P\{u) = P t ,(u) and P 2 {u) = P b ~(u). For 
the domain predominantly bound in well 2, we have P 2 (u) = P b + (u) and P\{u) = P b {u). 
For u 3> u c , we have P b ~ (u) ~ (see fig. pi a)), which implies that domain 1 then contains 
only the complexes R±Li, and domain 2 only the complexes R2L2. 

Since the total numbers Ni and N 2 of receptors 1 and 2 in the cell membrane are 
constant, we have 

Nx = [Ri]A + KxlR-RL^Ac (<j>P b + (u) + (1 - cj>)P b (u)) (24) 

N 2 = [R 2 ]A + K 2 [R 2 ][L 2 ]A C ((1 - <p)P b + {u) + (f>P b (u)) (25) 

Here, A is the total area of the cell, A c the contact area, and is the fraction of the 
contact area occupied by domain 1, which is predominantly bound in well 1. The first 
terms on the right-hand sides, [R\]A and [i^lA are the total numbers of unbound 
receptors. The concentrations [Ri] and [R 2 ] of unbound receptors within and outside 
of the contact zone are equal since we neglect the area occupied by bound receptor- 



ligand complexes within the contact zone (see also eq. (11)). The second terms on the 



right-hand sides of eqs. (24) and (25) are the numbers of bound receptors. In analogy 



to eq. (m), the number of bound receptors R± in the domain predominantly bound in 
well 1 is Ki[Ri][Li]A c <f)P^(u), and the number of bound receptors R 2 in this domain 
is K 2 [R 2 ][L 2 ]A c <pP^(u). The number of bound receptors R\ in the domain that is 
predominantly bound in well 2 is iiTi[i?i][Li]A c (l — 0)P 6 ~(«), and the number of bound 
receptors R2 in this domain is ir 2 [i?2][-^2]^.c(l — 4 ) )Pb~( u )- Below the critical potential 



P b (u) = Pb(u). The two equations (24) and (25) therefore 



depth u c , we have P b + (u) 
are independent from each other for u < u c , but dependent on each other for u > u c 
With the four independent, dimensionless parameters 



n\ 



NtKl 



2 

we 



A c k B T' 



n 2 



A c k B T 



and 



k 1 = K 1 [L 1 



I 



k 2 = K 2 [L 2 



the eqs. (|24|) and (]25|) can be rewritten as 

+ (1- 



m 



u 

u 



+ u(<t>P+ 



u 



A 

A 



)Pf(u)) 



n 2 = - r +u((l- cj>)P+(u) + <j)P b -{u)) 



(26) 

(27) 

(28) 
(29) 



since we have u = lfi[i?i][Li]^ e K/A; B T = K 2 [R2\[L2}ll e n/k B T (see eqs. Q to (g). 
From these two equations, one can determine u and <fi as functions of the independent 
parameters m, n 2 , k\ and k 2 . Domain coexistence in the cell contact zone occurs for 
u > u c and < (f> < 1. 

To obtain general relations for the critical point and the boundary lines of the two- 



phase region in the ki-k 2 plane, we first solve eq. (28) for k\ and eq. (29) for k 2 , which 
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Figure 6. (a) Monte Carlo data for the membrane fraction P& within the wells of a 
symmetric double- well potential with rescaled well depth u (see eq. (23 1). The data 
are from simulations with the rescaled width z wc = (/ W c/») V k/IzbT — 0.5 and rescaled 
separation Zba = {Iha./ a) ^k/IzbT — 2 of the two wells. Here, l wc and lb a are the width 
and separation of the wells, a is the linear size of the discrete membrane patches, and 
k is the effective rigidity of the membranes. Below the critical well depth u c w 0.053, 
the membrane is bound in both wells with the same fraction P&(u). Above the critical 
well depth u Cl the membrane is predominantly bound in one of the two wells. The 
membrane fraction bound in the dominant well is Pj. ,(u) (upper branch for u > u c ), 
and the membrane fraction bound in the other well is P^~(u) (lower branch for u > u c ). 
We have obtained the data from simulations with a square lattice of up to 200 x 200 
membrane patches and with up to 5 • 10 7 Monte Carlo steps per lattice site. Details of 
the Monte Carlo simulations are described in 18 21 . - (b) Exemplary phase diagram 
for the rescaled receptor numbers ri\ = 0.07 and ri2 — 0.09 (see eq. (26l) obtained from 



interpolation of the Monte Carlo data in (a) and insertion of the functions P^(u) and 



P b (u) in eqs. (33) and (34). Here, k\ and k^ are the rescaled concentrations of the two 



ligands in the supported membrane (see eq. (27)). The two-phase coexistence region is 
bounded by two lines along which the area fraction (f> of domain 1 in the contact zone 
= and ( 



is 



1 (see eqs. (33) and (34)) 
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leads to 



k 2 



a 



n 1 -u(0P+(u) + (l-0)p-(n)) 
u 

n 2 - u ((l-0)P+( u ) + 0p-( u )) 



(30) 
(31) 



At the critical point, we have P b + (u c ) = P b (u c ) = Pb(u c ). By inserting these relations 
into eqs. (30) and (31), we obtain a general expression for the location (ki, k 2 ) c of the 



critical point in the ki-k 2 plane: 



(h,k 2 



u r . 



U r 



n ± - u c P b (u c ) n 2 - u c P b {u c 



(32) 



For the Monte Carlo data of fig. |6[a) with u c ~ 0.053 and Pb(u c ) ~ 0.18, for example, we 
have {k 1 ,k 2 ) c ~ (0.053/ (n x - 0.0095), 0.053/(n 2 - 0.0095)). Since h x and k 2 have to be 
positive, domain coexistence can only occur if n\ and n 2 are both larger than u c Pb(u c ). 
The domain-coexistence region in the ki-k 2 plane is bounded by two lines with = 
and ( 



1. Inserting = 
{ki, ^2)0=0 



in the eqs. (30) and (|31|) leads to the parametric form 
u 



a 



m 



uP b {u) 



n 2 -uP b + {u 



for u > u r 



(33) 



for the = line. Similarly, inserting 
parametric form 

u 



1 in the eqs. (30) and (31) leads to the 



(* 



1; ^2 )</>=! 



■(/, 



m 



uP b + (u) 



n 2 - uP h (u) 



for u > u r 



(34) 



1 line. The domain-coexistence region in the phase diagram of fig. [6 



for the 4 

for example, follows from inserting the functions P b + {u) and P b {u) obtained from 

interpolation of the Monte Carlo data shown in fig. p[a) into eqs. (33) and (34). 



For u 3> u c , we have P b (u) w 0. The = line then is given by 



k\U\ 



4,=o n 2 - hn^ {kiUij 



since we have u 



k\ri\ for P b iu) w 0. Similarly, the = 1 line is given by 
k 2 n 2 



(35) 



(36) 



=1 m - k 2 n 2 P b (k 2 n 2 ) 
because of u ~ k 2 n 2 . For u 3> u c , the membranes are only bound via one of the wells. 
The membrane fraction P b bound in this well therefore can be approximated by the 
same expression P b {u) ~ u/{c\ + u) with c\ ~ 0.071 as in the case of an effective 
single- well adhesion potential (see eq. Q). 

The = line has a vertical asymptote in the k\-k 2 plane, since k 2 in eq. (35) 
diverges for 

k ini P b + {hm) (37) 



n 2 



because the denominator of the right-hand side of eq. (31) vanishes. With P l 
u/(ci + u), we obtain the location 

2ni V V n 2 



\u 



h 



(38) 
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for this vertical asymptote from eq. (37). Similarly, the = 1 line has a horizontal 



asymptote in the k\-k 2 plane since the the denominator of the right-hand side of eq. (36) 
vanishes for 

m = k 2 n 2 P^ (k 2 n 2 ) (39) 

With P^{u) ~ u/{c\ + u), we obtain the value 



^StW 1 ^ 1 



for the horizontal asymptote of the = 1 line from eq. (39) 



6. Discussion and Conclusions 

In this article, we have determined phase diagrams for cells adhering to supported 
membranes with anchored ligands. For supported membranes with short and long 
ligands L\ and L 2 that bind to the same cell receptor R, coexistence of domains of 
L\R and L 2 R complexes in the cell adhesion zone only occurs for equal effective binding 
strength Xi[Li] = K 2 [L 2 ] of the complexes where K\ and K 2 are the binding equilibrium 
constants at appropriate membrane separations. The domain coexistence thus occurs 
along a line in the [Li]-[L 2 ] plane, which ends at the critical point (see fig. |5](a)). We 
have assumed that the area of the supported membrane is much larger than the adhesion 
zone, which implies that the concentrations [Li] and [L 2 ] of unbound ligands do not 
change significantly during adhesion since the supported membrane constitutes a large 
'ligand reservoir'. Constant concentrations of unbound ligands imply constant chemical 
potentials /ii ~ /cBTln(a 2 [Li]) and fi 2 ~ kBT\n(a 2 [L 2 ]) of the ligands in our model (see 
eq. (32) in ref. (21|). A coexistence line as in the diagram of fig. [51(a)) is typical for 
phase diagrams in grand-canonical ensembles with constant chemical potentials. 

For supported membranes with two types of ligands L\ and L 2 that bind to different 
cell receptors R\ and R 2 , we obtain a qualitatively different phase diagram with a broad 
coexistence region (see fig. [5tb)). Domain coexistence in the cell adhesion zone occurs 
for i^i[i?i][Li] = K 2 [R 2 ] [L 2 ] . The broad coexistence region is a consequence of the fact 
that the concentrations of unbound receptors [Ri] and [R 2 ] depend on the numbers of 
bound receptors and, therefore, on the fractions of the cell adhesion zone occupied by the 
domains of R\L\ and R 2 L 2 complexes, since the total numbers N\ and N 2 of receptors 
in the cell membrane are constant. A broad coexistence region as in the diagram of 
fig. pVo)) is typical for phase diagrams in canonical ensembles with constant particle 
numbers. 

Milstein and coworkers (5] have observed domain coexistence for a narrow 
concentration ratio of short and long ligands that bind to the same cell receptor CD2, 
in agreement with our phase diagram in fig. [5ta). However, two differences between 
our phase diagram in fig. [51(a) and the phase diagram of Milstein and coworkers in 
fig. 7 of ref. [5] are: First, the coexistence line in the phase diagram of Milstein and 
coworkers seems to have a finite width. Such a finite width may result from slight 
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changes of the ligand concentrations upon binding, since several cells adhere to the 
same supported membrane in the experiments. Second, the coexistence line in the 
diagram of Milstein and coworkers ends in a region in which the cells do not adhere, 
while the coexistence line in the diagram of fig. [51(a) ends at a critical point. In this 
article, we have neglected repulsive interactions from, e.g., the cell glycocalyx. In our 
model, such repulsive interactions lead to an unbinding of the membranes at certain well 
depths Uf and C/| f of the effective adhesion potential shown in fig. B 20 . We obtain a 



phase diagram similar to the diagram of Milstein and coworkers if the well depths Uf 
and Iff at which the membranes unbind are larger than the critical potential depth Uf, 
which determines the location of the critical point in the diagram of fig. |5la). 

We find that thermal membrane shape fluctuations on nanometer scales play a 
central role during cell adhesion. Fluctuations on these scales have been recently 



reported for immune cells adhering to coated substrates 69 70 . In previous work, we 



have found that the fluctuations lead to a cooperative binding of receptors and ligands 



see fig. \3\j 21 , and to a critical point for the segregration of long and short receptor- 



ligand complexes 19 , 20 . Our phase diagrams in fig. 5 are therefore qualitatively 



different from phase diagrams calculated under neglection of shape fluctuations 43 
The binding cooperativity of receptors and ligands arises since a receptor-ligand complex 
locally constrains the membrane shape fluctuations and facilitates the binding of nearby 
complexes. The binding cooperativity is thus closely related to the fluctuation-induced 



attractive interactions between bound receptor-ligand complexes 17 18 29 , 71 72 , which 



result from a suppression of membrane-shape fluctuations, similar to the fluctuation- 



induced interactions of rigid membrane inclusions 73 76 . 

We have neglected here the line tension of the domain boundaries, which may 
suppress the formation of small domains in the cell adhesion zone. In classical nucleation 
theory, the line tension leads to a threshold size for stable domains. Experimental 
observations of stable microdomains in the adhesion zones of immune cells |9| [13|[65|[67] 
indicate that this threshold size is rather small. We will consider the line tension between 
domains of short and long receptor-ligand complexes in detail in a future article. 
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